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Magnetic response of interacting electrons in a zigzag carbon nanotube threaded by a magnetic 
flux is investigated within a Hartree-Fock mean field approach. Following the description of energy 
spectra for both non-interacting and interacting cases we analyze the behavior of persistent current 
in individual branches of a nanotube. Our present investigation leads to a possibility of getting a 
filling-dependent metal-insulator transition in a zigzag carbon nanotube. 
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I. INTRODUCTION 

The isolation of single layer graphene by Novoselov et 
at. [l| has initiated intense and diverse research on this 
system. Graphene, a single layer of carbon atoms tightly 
packed into a two-dimensional honey-comb lattice, has 
drawn attention of scientists in various disciplines due 
to its unconventional and fascinating electronic proper- 
ties arising particularly from the linear dispersion rela- 
tion around the Dirac points of the hexagonal Brillouin 
zone. These unique properties can be understood in 
terms of the Dirac Hamiltonian 0] since it actually de- 
scribes the physics of electrons near the Fermi level of 
the undoped material. The carriers in graphene effec- 
tively behave as massless relativistic particles within a 
low energy range close to Fermi energy and these mass- 
less Dirac Fermions Q evince various phenomena in this 
energy range. The bipartite character of the wonderful 
lattice structure of graphene strongly influences its intrin- 
sic properties and makes graphene a wonderful testbed 
not only for condensed matter theory, but also for quan- 
tum field theory and mathematical physics. Though lot 
of studies have been done both theoretically as well as 
experimentally to reveal electronic properties of this ex- 
otic system, yet complete knowledge about it is still lack- 
ing. This motivates us to address some interesting issues 
of electron transport in carbon nanotubcs where a nan- 
otube is formed by rolling up a graphite ribbon in the 
cylindrical form 

In this article we explore the behavior of persistent cur- 
rent in a finite sized graphite nanotube with zigzag edges 
within a nearest-neighbor tight-binding (TB) framework 
using a generalized Hartree-Fock (HF) approximation 0- 
0] ■ The phenomenon of persistent current in a mesoscopic 
ring, threaded by an Aharonov-Bohm (AB) flux <f>, has 
been studied theoretically more than two decades ago @- 
[TlT | . Later it has been justified through several nice exper- 
iments p^4l5j . It is a pure quantum mechanical effect 
and can persist without any dissipation in the equilib- 



rium case. In the last few years extensive studies on per- 
sistent current in carbon nanotubes have been performed 
and many interesting physical phenomena have been ex- 
plored |16l - [l8j . Persistent current in a carbon nanotube is 
highly sensitive to its radius, chirality, deformation, etc. 
Very recently it has also been observed experimentally 
that the Fermi energy of a carbon nanotube can be reg- 
ulated nicely by means of electron or hole doping, which 
can induce a dramatic change in persistent current fl8| . 
It is well established that in a conventional multi-channel 
mesoscopic cylinder electron transport strongly depends 
on the correlation among different channels as well as the 
shape of Fermi surface. Therefore we might expect some 
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FIG. 1: (Color online). Schematic view of a zigzag graphite 
nano-ribbon with N x and N y number of atomic sites along 
the x and y directions, respectively. 

interesting features of persistent current in a carbon nan- 
otube due to its unique electronic structure. 

The behavior of persistent current in zigzag nanotubes 
has been addressed theoretically by some groups fl6l - fl8j . 
There are many theoretical techniques available in the 
literature [L9l - |24| which generally investigate magnetic 
response of the entire system but the distribution of per- 
sistent current among different branches of the nanotube 
remains unaddressed though it is highly important to 
illuminate the magnetic response of the system with a 
deeper insight. To the best of our knowledge, no rigor- 
ous effort has been made so far to unravel the behavior 
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of persistent current in separate branches of a graphite 
nanotube. This is the main motivation behind this work. 

In what follows, we present the results. Section II is de- 
voted to present the model and generalized HF approach. 
Following the energy spectra for the non-interacting and 
interacting cases (Section III), in Section IV we estab- 
lish the second quantized form to evaluate persistent cur- 
rent in individual branches of a zigzag carbon nanotube. 
The energy-flux characteristics are described in Section 
V, while the behavior of persistent current in separate 
branches of a nanotube is illustrated in Section VI. Fi- 
nally, in Section VII we draw our conclusions. 



II. THE MODEL AND THE MEAN FIELD 
APPROACH 

We begin with a graphite nano-ribbon of zigzag edges 
as shown in Fig. [TJ where the filled magenta (large) 
and blue (small) circles correspond to two different sub- 
lattices, namely, A and B, respectively. N x and N y cor- 
respond to the number of atomic sites along the x and 
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FIG. 2: (Color online). A graphite nanotube threaded by a 
magnetic flux $. 

y directions, respectively. In order to elucidate magnetic 
response of a nanotube we roll up the graphite ribbon 
along x direction using periodic boundary condition and 
allow to pass a magnetic flux <j> (measured in unit of el- 
ementary flux quantum 0o = ch/e) along the axis of 
the tube as shown in Fig. [3J We describe our model 
quantum system by the nearest-neighbor TB framework 
which captures most of the essential properties of the 
tube nicely (25l - [28j . In presence of magnetic flux </>, the 
Hamiltonian of an interacting zigzag nanotube reads, 

H = t ^ (aJ ra>J j .6 m _i jn:(T e + al n n a .b m +i tn ,,j e 

m,n,cr 

+ a L«,<An,™+l^) + h.C. 

+ ( a L.n.t a ™."t a L.7U. a "',"4 



where, m and n are integers describing the co-ordinates of 
the lattice sites. The site indexing is schematically shown 
in Fig.[3]for better viewing, t is the nearest-neighbor hop- 
ping integral, n (b^ n ) is the creation operator for an 
electron of spin a (t,i) associated with A (B) type of 
sites at the position (m,n) and the corresponding anni- 
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FIG. 3: (Color online). Schematic view of different atomic 
sites with their co-ordinates. 

hilation operator is denoted by a mj „ (b m ,n)- The factor 
6 (= 2ircf>/N x ), the so-called Peierl's phase factor, is in- 
troduced into the above Hamiltonian to incorporate the 
effect of magnetic flux applied along the axis of the tube. 
U is the strength of on-site Hubbard interaction. 

Decoupling of interacting Hamiltonian: Using the gener- 
alized HF approach, we decouple the TB Hamiltonian 
into two different parts corresponding to two different 
values of a (j" and I). After decoupling, the Hamiltonian 
looks like, 

H MF = + H± + M o (2) 

where, 

= U 2^ {( n m,n,l) n m,n,t + ( n m+l,n,4.) n m+l,n,t) 
m,n 

+ { a LnA b rn-l,n^~ t0 + a^ nA b m+1 ^ t e W 

m, n 

+ a Ln,t 6 ™." + !-t + h ' C -) > ( 3 ) 

H\ = U ^ {( n m,n,-t) n m,n,i + ( n m+l,n,t) n m+l,n,i) 
m,n 

+ ( a L,n,i & m-l,n,4.e- ie + <4,„ 4 WmU ^ 

+ 4*,ra 6 ™,«+u + h - c -) (4) 

and 

+ ( n m+l,n,f) ( n m+l,n,i)) ■ (5) 

Here, n a and n b m n a are the number operators as- 
sociated with the A and B types of atoms, respectively, 
iff and H± are the Hamiltonians for up and down spin 
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electrons, respectively. Hq is a constant term which gives 
the energy shift. 

Self- consistent procedure: In order to get the energy 
eigenvalues of the interacting Hamiltonian we go through 
a self-consistent procedure considering initial guess val- 
ues of (n^ no .) and (ji^ nna ). With these initial values, 
the up and down spin Hamiltonians are diagonalized nu- 
merically and a new set of values of (n^ na ) and (n h m n a ) 
are calculated. These steps are repeated until a self- 
consistent solution is achieved. 

Finding the ground state energy: After getting the self- 
consistent solution we determine the ground state energy 
(Eq) at absolute zero temperature (T = OA") for a par- 
ticular filling by taking the sum of individual states upto 
the Fermi level (Ep) for both up and down spin electrons. 
The expression for ground state energy reads, 

E = j2 E i,t+J2 Ei >± +H ° ( 6 ) 

i i 

where, i runs over the states up to the Fermi level. -E^f's 
and Ei^s are the single particle energy eigenvalues ob- 
tained by diagonalizing the up and down spin Hamilto- 
nians iJ| and H±, respectively. 



III. ENERGY SPECTRUM 

To make this present communication a self contained 
study let us first start with the energy band structure of 
a finite width zigzag nano-ribbon. 

Non-interacting case: To establish the energy dispersion 
relation of a zigzag nano-ribbon we find an effective dif- 
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FIG. 4: (Color online). Unit cell configuration of a zigzag 
nano-ribbon. 

ferencc equation analogous to the case of an infinite one- 
dimensional chain. This can be done by proper choice of 
a unit cell from the nano-ribbon. The schematic view of 
a unit cell configuration with N y pairs of B-A atoms in a 
zigzag nano-ribbon is shown in Fig. 2) With this arrange- 
ment, the effective difference equation of the nano-ribbon 



gets the form, 

(El - £ a )^ a = Tlp j+ X, a + 7tyi-l,«r (?) 

where, 



V ^jN y A,a I 

£ and T are the site-energy and nearest-neighbor hopping 
matrices of the unit cell, respectively. I is a (2N y x 2N y ) 
identity matrix. Since in the nano-ribbon translational 
invariance exists along the rr-direction, we can write ipj^ 
in terms of the Bloch waves and then Eq. [7] takes the 
form, 

(EI - £ a ) = Te lk * K + 7-t e -tt-A (9) 

where, A = \/3a is the horizontal separation between two 
filled magenta or blue circles situated at two successive 
unit cells, a is the length of each side of a hexagonal 
benzene like ring. Solving Eq.[9]we get the desired energy 
dispersion relation (E vs. k x ) of the ribbon. 

As illustrative example, in Fig. [S] we show the variation 
of energy levels (green curves) as a function of wave vec- 
tor k x for a finite width zigzag nano-ribbon considering 
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FIG. 5: (Color online). Energy levels (green curve) as func- 
tion of k x for a finite width zigzag nano-ribbon considering 
N y — 4. The discrete eigenvalues (filled black circles) of a 
nanotube with N x = 12 and N y = 4, in the absence of AB 
flux <j>, are superimposed. Here we set U = 0. 

N y = 4. Quite interestingly we observe that at E = 0, 
partly flat bands appear in the spectrum which make 
the system unique. The electronic states corresponding 
to those almost flat bands are characterized by strongly 
localized states near the zigzag edges of the tube. The 
existence of these edge states have also been reported 
earlier by some other groups 0, HH, H3| • 

Following the energy band structure of a finite width 
nano-ribbon now we focus on the variation of energy 
levels of a nanotube. For a nanotube k x also be- 
comes quantized where the quantized values are obtained 
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by applying periodic boundary condition along the x- 
direction [31| . The quantized wave numbers are ex- 
pressed from the relation k x = 4nn x /N x A, where n x is 
an integer lies within the range: —N x /A < n x < N x /A. 
Plugging the quantized values of k x in Eq. |9]we can easily 
determine the eigenvalues of a finite sized nanotube. As 
representative example, in Fig. [5] we show the variation 
of discrete energy eigenvalues (filled black circles) for a 
zigzag nanotube considering N x = 12 and N y — 4, in the 
absence of AB flux <fi passing through the tube. For this 
nanotube k x gets six quantized values (— 7r/A, — 27r/3A, 
— 7r/3A, 0, 7r/3A and 27r/3A), and therefore, total 48 en- 
ergy values are obtained since N y is set at 4. 

Interacting case: In the presence of e-e interaction energy 
levels get modified significantly depending on the filling 
of electrons. The results calculated for a particular value 



v<t {=x a ) as, 




(10) 



where, x a is the displacement operator. The velocity 
operator is computed from the expression, 



in 



(11) 



Using this relation we can write the velocity operator of 
an electron with spin a in a zigzag channel n (say) in the 




FIG. 7: (Color online). Energy-flux characteristics of a half- 
filled zigzag nanotube with N x = 10 and N y = 7. (a) U = 
and (b) U = 1.5. 



FIG. 6: (Color online). Energy levels as function of k x for 
a finite width zigzag nano-ribbon considering N y = 4 and 
U = 1.4, where (a) and (b) correspond to the one-third- and 
half-filled cases, respectively. 

of U are presented in Fig. |6] where we set N y = 4. In 
the half-filled band case, a gap opens up at the Fermi 
energ y 1321 which is consistent with the DFT calcula- 
tions [331 ] and the gap increases with the value of U. 
A careful investigation also predicts that the full energy 
band gets shifted by the factor U/2. 



IV. SECOND QUANTIZED FORM OF 
PERSISTENT CURRENT 

In order to evaluate persistent current in individual 
zigzag paths of a nanotube, threaded b y a n AB flux </>, 
we use second quantized approach [34l |35|. This is an 
elegant and nice way of studying the response in separate 
branches of any quantum network. 

We start with the basic equation of current operator / CT 
corresponding to spin a in terms of the velocity operator 



form, 



4- a) b i e' w ) 



(12) 



Therefore, for a particular eigenstate \ip p>a -) persistent 
current in n-th channel becomes, 



(13) 



where, the eigenstate \ip p , a ) looks like, 



+ P P m+l,n,*\ m + I,": ') 



(14) 



Here, |m, n, a)'s are the Wannier states and aP m n a and 
@mna s arc ^hc corresponding coefficients. Simplifying 
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Eq. [T3l we get the final relation of persistent charge cur- Using the same prescription we can also evaluate per- 

rent for n-th zigzag channel as, sistent current in individual armchair paths (along y di- 

rection) of the nanotube. The final expression of it gets 



_ le -t \ " (op * p -ie the form 



n p * rp _ a. 

m,n,<7^m+l,7i,(7 ' 



P * rvP e i9 



m— l,n,cr m,n,(r 
v p * B p e~ ie ) ClM 

*777, 71,(7^771 — 1,71,(7 / ' \ J 
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J m-1, 



it 



2hN„ 



N v N y 

E(ry p * B p p~ ie - B p * rrrr p r W ) 4- ( B p * n p e~ ie - n/ p * 

V 777,77, CrMm — 1,71,(7 "777 — 1,77' 777,71,(7^ ) 1 / j V^TTl ,77 ,(7^771 — 1 ,77,(7 777—1,77,(7 

77=1,3,... 



77=2,4,... 

jV„-l 



X /j p p i61 4- B p * rt p - rv p * B p ) 4- ( B p * rv p - rt p * B p ) 

" ym.n^a^ ' t^rn — 1,71,(7^771 — 1,71 — 1,(7 771—1,77—1,(7^771 — 1,71,(7/ ^ / j \r VTi ,71+ 1 ,<T ^777 ,71, (7 ^ 771,77 ,(7 M777 ,71+ 1 , (7 / 



71=2.4, 



(16) 



where, an armchair channel (to — l,m) is constructed by 
(to — l)-th and TO-th lines according to our indexing. 

At absolute zero temperature (T = OK), net persis- 
tent current driven by electrons of spin a in a particular 




FIG. 8: (Color online). Few energy levels of Fig. [7] are re- 
plotted for a narrow energy range across the band centres for 
better viewing of the variation of energy levels with flux 4>, 
where (a) and (b) correspond to the identical meaning as in 

Fig.m 

channel n for a nanotube described with Fermi energy 
Ep can be determined by taking the sum of individual 
contributions from the lowest energy eigenstates upto the 
Fermi level. Hence we get, 



(17) 



Summing l n , a over all possible channels n, and a we get 
total persistent current in the nanotube which is mathe- 
matically expressed as, 



71,(7 



(18) 



The persistent current can also be determined in some 
other ways as available in literature. Probably the sim- 




FIG. 9: (Color online). Energy gap {AE) as a function of 
on-site Hubbard interaction strength U for a zigzag nanotube 
with N x — 10 and N y — 7 in the half-filled band case when cj> 
is set at </>o/2. 

plest way of determining persistent current is the case 
where first order derivative of ground state energy with 
respect to AB flux <fi is taken into account. Therefore, 
we can write, 



dE Q (<f>) 



(19) 



where, Eo(<fi) is the total ground state energy for a par- 
ticular electron filling. But, in our present scheme, the 
so-called second quantized approach, there are some ad- 
vantages compared to other available procedures. Firstly, 
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we can easily measure current in any branch of a compli- 
cated network. Secondly, the determination of individual 
responses in separate branches helps us to elucidate the 
actual mechanism of electron transport in a more trans- 
parent way. 

In the present work we examine all the essential fea- 
tures of persistent current at absolute zero temperature 
and use the units where c = h = e = 1. Throughout our 
numerical calculations we set t = — 1 and measure all the 
physical quantities in unit of t. 



V. ENERGY-FLUX CHARACTERISTICS 

In Fig. [Jj we show the variation of energy levels as 
a function of flux <p for a zigzag nanotubc considering 
N x = 10 and N y = 7 both for the (a) non-interacting 
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FIG. 10: (Color online). Ground state energy level as a func- 
tion of </j for a zigzag nanotube in the half-filled band case 
considering N x = 20 and N y = 8. (a), (b), (c) and (d) corre- 
spond to U = 0, 0.5, 1 and 1.5, respectively. 

and (b) interacting cases. For U = 0, we compute the 
energy levels simply by diagonalizing the non-interacting 
Hamiltonian and the nature of the energy spectrum be- 
comes independent of the total number of electrons N e in 
the system. On the other hand, for the non-interacting 
case (U 7^ 0) we fist decouple the interacting Hamiltonian 
(Eq. [1} for a particular filling, in the mean field scheme, 
into two non-interacting Hamiltonians (for up and down 
spin electrons) and then diagonalize the Hamiltonian for 
up (down) spin electrons. For identical filling factor of 
up and down spin electrons the energy levels are exactly 
similar both for and (see Fig. [7[b)), and there- 
fore, one energy spectrum cannot be separated from the 
other. Since in our case we set N x = 10 and N y = 7, 
we get total 70 independent energy levels and due to 
their overlaps individual energy levels are not clearly dis- 
tinguished from the spectra given in Fig. [JJ To have a 
better look in Fig.|5]we re-plot few energy levels of Fig. [JJ 
collecting them from a narrow energy range across the 



band centres, where (a) and (b) correspond to the iden- 
tical meaning as in Fig. [7J From the spectra we see that 
all the energy levels vary periodically with providing 
4>o (= 1 in our chosen unit) flux-quantum periodicity. 
At half- integer or integer multiples of <pg, energy levels 
have either a maximum or a minimum (see Fig. [5]) , and 
accordingly, at these points persistent current becomes 
zero which is quite obvious since the current is obtained 
by taking the first order derivative of the eigenenergy 
with respect to flux <f> (Eq. \T§§. Both the energy spec- 
tra take a complicated look as there are many crossings 
among the energy levels particularly in the regions away 
from E = (Fig. [8]). At E = 0, the energy levels become 
almost flat for a wide range of <f>, and, near <p = ±</5o/2 
they vary slowly with </5 as shown in Fig. HJa). These al- 
most flat energy levels support a very little contribution 
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FIG. 11: (Color online). Current-flux characteristics of a 
zigzag nanotube with N x — 14 and N y — 6, where (a), (b), (c) 
and (d) correspond to N e — 10, 20, 30 and 82, respectively. 
U is fixed at 1. 

to the persistent current as clearly followed from Eq. Il~9l 
While, the other energy levels with larger slopes provide 
large persistent current. This peculiar nature of the en- 
ergy levels invokes the current amplitude to become filling 
dependent and we elaborate it in the following section. In 
Fig. [Jjb) we display the variation of energy levels with 
<fi for a zigzag nanotube with the same parameter values 
declared above in the presence of Hubbard interaction. 
Here we choose U = 1.5. Both for the up and down 
spin Hamiltonians the eigenvalues are exactly identical 
and they overlap with each other. Few energy levels of 
Fig. EJb) are also re-plotted in Fig. [SJb) for better view- 
ing. The electronic correlation leads to an energy gap 
at the band centre and the gap increases with U. It is 
illustrated in Fig. [5J This energy gap is consistent with 
the energy gap obtained in the E-k x diagram (Fig.[5Jb)). 

The variation of ground state energy level of a carbon 
nanotube with zigzag edges as a function of magnetic 
flux <fi is depicted in Fig. [TU]in the half-filled band case, 
where (a), (b), (c) and (d) correspond to the four differ- 
ent values of electronic correlation strength U = 0, 0.5, 1 
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FIG. 12: (Color online). Persistent current in individual zigzag paths as a function of 4> for a half-filled zigzag nanotube 
(N x — 20 and N y = 7) with U = 1.2, where, (a)-(g) correspond to lst-7th zigzag channels of the tube, respectively. The net 
current corresponding to both up and down spin electrons are displayed in (h) while in (i) total persistent current is shown. 



and 1.5, respectively. The energy levels evince one flux- 
quantum periodicity, as expected, and their energies get 
increase with U. It can be explained as follows. In pres- 
ence of U, each site is occupied by an electron with either 
up or down spin in the half-filled band case. This is the 
ground state energy configuration. Now, if we add an 
electron further, the probability of getting two opposite 
spin electron in a single site becomes finite which gives 
higher energy due to the repulsion in presence of U. 

VI. CURRENT-FLUX CHARACTERISTICS 

Now we focus our attention on the behavior of persis- 
tent current in a zigzag nanotube. 

First, we illustrate the dependence of persistent cur- 
rent amplitude on the electron filling. To ensure it in 
Fig. [11] we present the current-flux characteristics of a 
zigzag nanotube considering N x = 14 and N y = 6, where 
four different figures correspond to the four different cases 
of electron filling. The Hubbard interaction strength is 
set at 1. It is observed that when the number of elec- 
trons is much smaller than half-filling, persistent cur- 



rent exhibits multiple kinks at different values of 0, as- 
sociated with the multiple crossings of energy levels, as 
shown in Figs. [Tlla'l-fcV In these three cases the num- 
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FIG. 13: (Color online). Total persistent current obtained in 
a traditional derivative approach (Eq. I19p as a function of <j> 
for the same parameter values mentioned in Fig. 1121 

ber of electrons are 10, 20 and 30, respectively. This is 
quite analogous to the nature of persistent current ob- 
served in conventional multi-channel mcsoscopic cylin- 
ders. The behavior of persistent current gets significantly 
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modified when the nanotube becomes half-filled or nearly 
half- filled. For example, see Fig. [TTT dV Here we choose 
N e = 82 i.e., the nanotube is very near to the half- filled 
band case. In such a case all the kinks disappear and 
current varies almost continuously, analogous to the be- 
havior of persistent current observed in traditional single- 
channel mesoscopic rings. For the cases when the nan- 
otube is far away from half-filling, current amplitudes are 
quite comparable to each other (see Figs.HlTa)-(c)). On 
the other hand, when the tube is nearly half-filled current 
amplitude remarkably gets suppressed. It is illustrated 
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FIG. 14: (Color online). Persistent current in an armchair 
path as a function of cj> for a zigzag nanotube with the same 
parameter values as mentioned in Fig. 1121 

in Fig. Illf d) . This enormous reduction of current ampli- 
tude can be visualized from the E-<p characteristics given 
in Fig. [7] At half-filling or very close to half-filling, the 
top most filled energy level lies in the nearly flat region 
i.e., around E = (sec Fig. [2a)) and it contributes a 
little to the current. Moreover, when U ^ 0, there is 
gap in the midband region. Now, for a particular filling 
the net persistent current is obtained by taking the sum 
of individual contributions from the lowest filled energy 
levels, and, in this process only the contribution which 
comes from the highest occupied energy level survives fi- 
nally and the rest disappear due to their mutual cancel- 
lations. It leads to the enormous reduction of persistent 
current amplitude in the half-filled or nearly half-filled 
case. This feature is independent of the size of the nan- 
otube. In all these cases, persistent current varies period- 
ically with flux <p, exhibiting 0o flux-quantum periodicity. 
From these current-flux spectra we can emphasize that, 
the current amplitude in a zigzag nanotube is highly sen- 
sitive to the electron filling and this phenomenon can be 
utilized in designing a high conducting to a low conduct- 
ing switching operation and vice versa. 

Finally, we concentrate on the behavior of persistent 
current in separate branches of a zigzag carbon nanotube. 
As illustrative examples, in Fig. 1121 we show the variation 
of persistent current in individual zigzag paths as a func- 
tion of flux (j) for the half-filled case considering N x = 20 
and N y = 7, where (a)-(g) correspond to 1st- 7th zigzag 
channels of the tube, respectively. The Hubbard correla- 
tion strength U is set equal to 1.2. In each of these figures 
we display currents carried by both up and down spin 



electrons together. They are exactly superposed with 
each other. All these currents exhibit 0o flux-quantum 
periodicity and their magnitudes are quite comparable to 
each other. Interestingly we see that Ii^ is exactly iden- 
tical to Irf, and, similarly for the (l2f, let) an d C^t) ^5f ) 
pairs. /4,f, the current in the middle channel, becomes 
the isolated one since we have chosen N v = 7. This is 
true for any zigzag nanotube with odd N y . For a tube 
with even N y , currents are pairwise identical. Summing 
up the individual currents in seven zigzag channels we get 
net persistent current carried by up and down spin elec- 
trons in the nanotube which is presented in Fig. ll2f h) and 
the total current is displayed in Fig. PT2T i) which exactly 
matches with the total current derived from the conven- 
tional method where first order derivative of the ground 
state energy is taken into account, as shown in Fig. 1131 
It emphasizes that the net contribution of persistent cur- 
rent in a zigzag carbon nanotube comes only from the 
individual zigzag channels, not from the armchair paths. 
To justify it in Fig. Q3]we present the variation of persis- 
tent current in an armchair path as a function of <j) for 
a half-filled zigzag nanotube considering N x = 20 and 
N y = 7, which clearly shows zero current for the entire 
range of <fi. 



VII. SUMMARY 

To conclude, in the present work we investigate in de- 
tail the magnetic response of a zigzag nanotube, threaded 
by a magnetic flux 0, using a generalized Hartree-Fock 
mean field approach. The model is described by a sim- 
ple tight-binding framework. Following the E-k x spectra 
of both the non-interacting and interacting cases of a 
finite width nanoribbon, we present the results of a nan- 
otube with zigzag edges. Energy levels get modified sig- 
nificantly in the presence of Hubbard interaction and the 
nature of the energy spectrum strongly depends on the 
electron filling. At the half-filled bad case, a gap opens 
up at the Fermi energy which is consistent with the DFT 
calculations. Next, we establish the second quantized 
form to evaluate persistent current in individual paths of 
a zigzag carbon nanotube. From the current-flux char- 
acteristics we can emphasize that the current amplitude 
in the zigzag nanotube is highly sensitive to the electron 
filling and this phenomenon can be utilized in designing 
a high conducting to a low conducting switching device 
and vice versa. 
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